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摘要 : 鼓 泡 型 板式 换 热 器 是 一 种 新 型 热 交 换 器 。 目 前 其 板 型 几何 参数 对 板 片 通道 内 流动 与 传 热 的 影响 规律 还 不 清楚 。 本 文 基于 
多 场 协同 原理 数值 研究 了 板 型 主要 几何 参数 对 流动 传 热 的 影响 规律 ， 发 现 当 单 元 体 长 度 增 大 或 宽度 减 小 时 ， 通 道 摩擦 阻力 因子 
f 减 小 ， 综合 性 能 评价 参数 PEC 值 增 大 ; 当 单 元 体 相 对 主流 流动 方向 倾斜 角 为 90" 时 ， f 值 最 小 ， PEC 值 最 大 。 本 文 的 研究 结 
果 对 鼓 泡 板 片 传 热 元 件 的 几何 结构 优化 及 板 型 优化 设计 具有 指导 意义 。 
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The effect of the structure of bubble type plate element on heat transfer 
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Abstract: Bubble plate heat exchanger is a new type heat exchanger. The influence of the plate geometric parameters 
on the flow and heat transfer is unclear so far. Based on the principle of multi-field synergy, the influence of the plate 
geometric parameters on the flow and heat transfer is studied numerically. It is find that when the unit length 
increases or the width decreases, the friction factor f decreases, performance evaluation parameters PEC increases; 
when the inclination angle of the unit body relative to the main flow direction is 90 °, the factor f obtain the 
minimum value, PEC obtain the maximum value. The results can be applied to the optimization of the bubble plate 
heat transfer element and the plate design. 
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0 引言 


板式 热 交 换 器 具有 传 热 系数 高 、 阻 力 系数 小 、 的 协同 关系 来 评价 传 热 元 件 的 性 能 ， 这 些 理 i 
结构 紧 竣 、 运 行 费 用 低 等 特点 ， 被 广泛 运用 于 热能 ” ”为 “多 场 协 同 ” 原 理 。 
动力 、 石 油 化 工 等 领域 。 鼓 泡 型 板式 换 热 器 是 一 种 本 文 数值 模拟 了 鼓 泡 型 板 片 元 件 组 成 通道 内 
新 型 板式 换 热 器 ， 采 用 无 源 强化 传 热 技术 来 改善 通 的 流 场 ， 基 于 “多 场 协同 ”原理 对 鼓 泡 板 片 传 热 元 
道内 的 流动 与 传 热 性 能 ， 其 主要 优点 是 流动 分 配 均 。 ” 件 的 几何 参数 对 流动 传 热 的 影响 进行 了 评价 ， 得 出 
匀 ， 摩 擦 压 降 小 ， 板 片 单 位 面积 上 接触 点 少 ， 有 利 几何 参数 对 流动 传 热 的 影响 规律 ,得 到 了 优化 参数 。 
于 减少 结 垢 与 沉积 中 。 

文献 23 提出 了 “ 场 协同 ”原理 ， 揭 示 了 减 注 1 数值 模拟 方法 
热 边界 层 厚度 与 增加 流体 扰动 之 间 的 内 在 关系 ， 认 
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为 减 小 主流 速度 矢量 与 温度 梯度 矢量 之 间 的 协同 角 个 1 物理 模型 与 数学 描述 

可 以 强化 对 流 换 热 。 文献 7 建立 了 速度 场 与 主流 方 鼓 泡 板 片 主 要 由 单元 体 长 度 /宽度 b、 高 度 及 
向 速度 梯度 场 之 间 的 协同 关系 ， 发 现 流体 沿 主流 方 相对 主流 方向 的 倾斜 角 01 以 及 夹 角 久 五 个 基 
向 速度 矢量 与 温度 梯度 矢量 间 的 协同 关系 可 以 表征 ” 本 几何 参数 确定 ， 板 片 及 单元 体 结 构 如 图 1 所 示 。 


传 热 元 件 的 综合 性 能 。 上 述 学 者 通过 寻找 各 矢量 间 计算 域 如 图 2 所 示 ， 在 横 纵 方向 上 各 包含 4 个 单元 
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体 。 其中, h 为 2.8 mm，0 为 90” ，7! 取 值 分 别 为 
20、22、24、26、28、30(mm), b 取 值 分 别 为 9、 
10、11、12、13、14(mm) 及 91 的 取 值 分 别 为 0°、 


15°、30°、45°、60°、75°、90°, 板 长 工 为 169.48 mm， 
板 宽 历 为 149.48 mm。 


图 1 鼓 泡 型 板式 热 交 换 器 及 单元 体 结构 
Fig.1 Bubble type plate heat exchanger and its structure 
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图 2 鼓 泡 板 片 通道 计算 域 


Fig.2 Bubble plate channel computational domain 


鼓 泡 型 板 片 通道 类 似 于 长 方形 罕 通 道 ， 通 道内 
单 相 流动 临界 雷诺 数 小 于 圆 管内 流动 临界 雷诺 数 。 
一 般 ， 训 泡 板 片 通道 中 的 流动 处 于 低 雷 诡 数 清流 状 
态 &29， 本 文 计算 采用 的 雷诺 数 大 约 为 2200。SST 
淇 流 模型 综合 了 k-e 和 k-w 两 种 模型 的 特点 ， 近 壁 
面 区 使 用 kw 模型 预测 剪 切 流动 ， 主 流 区 使 用 ke 
模型 预测 汕 流 流动 4， 既 考 虑 壁面 附近 汕 流 前 切 应 
力 传输 ， 又 可 预测 流动 开始 和 负 压 力 梯度 下 流体 分 
离 ， 不 会 对 涡流 粘度 造成 过 度 预 测 ， 适 用 于 鼓 泡 型 
板 通道 内 的 流 场 计 算 。 


1. 2 边 宽 条 件 与 网 格 验证 

采用 ICEM 进行 非 结 构 化 四 面体 网 格 划 分 ， 壁 
面 用 三 棱柱 网 格 加 密 , 网 格 验证 采用 的 几何 尺寸 为 : 
三 20 mm，b=9 mm，01=0;=90°。 工 质 为 水 ， 入 口 为 
速度 入 口 边界 , 入口 温度 293K, 采用 压力 出 口 边界 ; 
壁面 为 无 滑 移 边 界 ， 板 壁面 等 壁 温 7,=333K， 通 道 
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两 侧 绝 热 。 采用 CFX 软件 , 网 格 独 立 性 验证 后 采 月 
了 6597193 个 网 格 ， 对 应 的 网 格 尺寸 为 0.72 mm， 
正式 计算 中 网 格 尺 寸 均 为 此 值 。 


1. 3 数据 处 理 方法 


摩擦 阻力 系数 采用 范 宁 摩擦 阻力 系数 : 
Ap 站. 
f= 5p Le (1) 
式 中 ，Zpr 为 通道 进出 口 间 的 流动 摩擦 压 降 ， 
DD。 表示 通道 当量 直径 日 De=47。 

壁面 局 部 Nuwan 数 及 壁面 平均 Nun 数 为 : 


hD 
Nuva = 2) 
| Ni UvandA 
Nu = G) 
Aan 


式 中 ，h 表 示 壁 面 局 部 对 流 换 热 系数 ，4 表示 
水 的 导热 系数 。 

通常 传 热 元 件 换 热 增强 ， 流 动 阻力 也 会 增 
大 。 这 里 使 用 如 下 综合 性 能 评价 参数 ， 
Re @) 
(7 万 
式 中 ，N 和 态 分 别 代表 光板 的 努 塞 尔 数 和 麻 
擦 阻力 因子 。 


1.4 多 场 协 同 原理 

基于 “多 场 协 同 ” 原 理 *"， 采 用 速度 场 矢 
量 UU 与 沿 主 流 方向 分 速度 梯度 场 撩 量 Vu 的 协 
同 角 a 评价 不 同 几 何 参 数 下 的 摩擦 阻力 特性 : 
L7 .VE 
TO (5) 
wa 

采用 速度 场 矢量 U 与 温度 梯度 场 矢量 VY7 
的 协同 角 8 评价 不 同 几 何 参数 下 的 传 热 特 性 ; 


U:.VvT 
0 = arccos 加 wz (0) 


采用 主流 方向 分 速度 梯度 场 天 量 Vu 与 温 
度 梯 度 场 矢量 VT 的 协同 角 评价 不 同 几何 参 
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CQ 三 G1CCOS 


数 下 的 综合 特性 : 
‘Vu 
Y= arccos v7 | [va (7) 
整 场 平均 协同 角 取 积分 平均 中 值 角 1: 
SU:vudr 
CQ 三 G1CCOS (8) 


DUl-lvular 


DO =arccos ed (9) 
> oflv7lar 
> V7.VudF 
NJ 三 QG1CCOS (10) 
v7|-|vular 


2 结果 与 讨论 


2.1 单元 体 长 度 /对 整 场 平均 协同 角 的 影响 

固定 单元 体 宽度 为 9 mm， 相 对 主流 方向 的 倾 
斜 角 01 以 及 结构 夹 角 0 均 为 900", 图 3 给 出 了 平均 
协同 角 随 单元 体 长 度 1 的 变化 规律 。 随 着 1 的 增 大 ， 
速度 场 U 与 分 速度 梯度 场 Vu 的 平均 协同 角 om 增 
大 ，am<90”， 二 者 的 协同 性 变 弱 ， 摩 擦 阻力 减 小 ， 
了 值 减 小 ; 速度 场 UU 与 温度 梯度 场 VT 的 整体 协同 角 
bm 减 小 ，Bm<90”， 二 者 的 协同 性 变 强 ， 传 热 效 果 
增强 ， 平 均 Nun 数 增 大 ;温度 梯度 场 VT 与 分 速度 


(@) 


梯度 场 Vu 的 平均 协同 角 知 增 大 ， 且 yw>90”， 二 
者 的 协同 性 变 强 ， 综 合 性 能 评价 参数 PEC 值 增 大 。 
2. 2 单元 体 宽度 b 对 整 场 平均 协同 角 的 影响 

固定 单元 体 长 度 1 为 20mm， 相 对 主流 方向 的 
倾斜 角 0 以 及 结构 夹 角 0, 均 为 90%, 图 4 给 出 了 平 


均 协同 角 随 单元 体 宽度 b 的 变化 规律 。 随 着 5 的 增 
大 ， 速 度 场 UU 与 分 速度 梯度 场 Vu 的 平均 协同 角 wm 


减 小 ，am<90”， 二 者 的 整体 协同 性 变 强 ， 摩 擦 阻 
力 增 大 , f 值 增 大 ;， 随 着 4b 的 增 大 ， 速 度 场 U 与 温 


度 梯度 场 VT 的 整体 平均 协同 角 pi 减 小 ， 
Bm<90”， 二 者 的 整体 协同 性 变 强 ， 传 热 效 果 增 强 ， 
平均 Nun, 数 增 大 ; 随 着 的 增 大 , 温度 梯度 场 VT 与 
分 速度 梯度 场 Vu 的 平均 协同 角 知 一直 减 小 ， 且 
ym>90”, 二 者 的 整体 协同 性 变 车， 综合 性 能 评价 参 
数 PEC 值 会 减 小 。 
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图 3 三 种 平均 协同 角 随 7 的 变化 规律 
Fig.3 Three average synergy angles variation with / 
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Fig.4 Three average synergy angles variation with b 


2. 3 单元 体 倾斜 角 01 对 整 场 平均 协同 角 的 影响 

5 所 示 为 91 变化 时 三 种 平均 协同 角 的 变化 规 
律 ， 其 中 三 20 mm,， b=9 mm，0, 为 90?。 图 5 (a)， 
当 91=30” 时 ， 通 道内 速度 场 U 与 主流 方向 分 速度 
梯度 场 Vu 的 平均 协同 角 wm 最 小 ， 且 am<90”， 二 


区 


者 的 整体 协同 性 最 强 ，f 值 最 大 ; 而 当 01=90” 时 ， 
am 最 大 ， 对 应 f 值 最 小 。 由 图 5 (b)， 随 0 的 变化 ， 
平均 协同 角 会 先 减 小 后 增 大 ， 且 当 0=45” 时 ， 
bun 最 小 ， 壁 面 平均 Nun, 数 最 大 。 由 图 5 (c)， 随 91 
的 增 大 , 平均 协同 角 知 先 减 小 后 增 大 , 且 yw>90”， 
当 0=90” 时， 温度 梯度 场 V7 与 分 速度 梯度 场 Vu 


的 


协 同性 最 强 ， ym 最 大 ， 综合 性 能 参数 PEC 值 最 大 。 
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